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We present a theory of nuclear magnetic resonance (NMR) shielding tensors for electronic states 
with arbitrary degeneracy. The shieldings are here expressed in terms of generalized Zeeman (g^) 
and hyperfine (A^) tensors, of all ranks k allowed by the size of degeneracy. Contrary to recent 
proposals [T. O. Pennanen and J. Vaara, Phys. Rev. Lett. 100, 133002 (2008)], our theory is valid 
in the strong spin-orbit coupling limit. Ab initio calculations for the 4-fold degenerate Pg ground 
state of lanthanide-doped fiuorite crystals CaF2:Ln (Ln = Pr 2+ , Nd 3+ , Sm 3+ , and Dy 3+ ) show that 
previously neglected contributions can account for more than 50% of the paramagnetic shift. 



Introduction. Paramagnetic nuclear magnetic reso- 
nance (pNMR) spectroscopy is a fundamental tool for 
probing static and dynamic local magnetic properties of 
materials [T] and metallo-proteins [2]. pNMR plays a 
central role in the elucidation of the quantum dynamics 
of single-molecule magnets and antiferromagnetic spin 
rings [3J. Moreover, it is an increasingly central tech- 
nique for probing strong spin-orbit coupled electronic 
states in a dissipative environment, as shown by a re- 
cent study of quantum tunneling processes [I], silent in 
ac-susceptibility experiments, in lanthanidc-based (Dy 3+ 
and Tb 3+ ) molecular nanomagnets. 

Despite the central role of pNMR in the development 
of new magnetic materials, only very recently ab initio 
approaches have been developed for the calculation of 
fundamental pNMR observables, such as nuclear shield- 
ing tensors [5J. In particular, by generalizing the work 
of Moon and Patchkovskii [5] on Kramers doublets, the 
paper by Pennanen and Vaara [7] stands out as the first 
comprehensive formulation of NMR shielding tensors in 
terms of molecular response properties, like the second- 
rank (first-rank in spin) EPR g and hyperfine A ten- 
sors, routinely computed via ab initio methods [5HTHTT)]. 
The work by Pennanen and Vaara [7] aims at proposing 
an NMR-theory for arbitrary electronic ground state de- 
generacy. However, this task is only accomplished in [7] 
for the weak spin-orbit coupling limit, and to date no 
theory fully accounts for truly general (spin and orbital) 
electronic degeneracy, as that arising in strong spin-orbit 
coupled lanthanide nanomagnets. 

In this Letter we propose a general theory for the ab 
initio calculation of pNMR chemical shifts as function 
of molecular tensors of all ranks consistent with the size 
of degeneracy, which applies to arbitrarily strong spin- 
orbit coupling (SOC) limit, and any degree of electronic 
degeneracy. 

Theory. In the linear response regime, the electron- 
induced magnetic field experienced by a nucleus I 
in a molecule immersed in a magnetic field B, is pro- 
portional to the field itself, so that B^ nd = — B • er 1 . 
This leads to an electron-mediated interaction energy be- 
tween external field and the magnetic dipole moment fi 1 
of nucleus / which is bilinear in the field and the nuclear 



dipole, given by W B ' ,X = — Bj r nd • fj, 1 . The proportional- 
ity tensor a 1 is the NMR shielding tensor, its trace being 
related to the chemical shift of nucleus /. 

In a system with a degenerate ground state, the en- 
ergy VP B,/X must be averaged among the thermally ac- 
cessible states resulting from the splitting of the de- 
generacy under the experimental conditions, leading to 
(W B >»') = -(Bf nd ) ■ y. 1 = B a {a^}^, where Greek in- 
dices correspond here to Cartesian tensor components, 
and a sum is implied over repeated indices. Then the 
field-independent shielding tensor is defined as 
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A strategy to evaluate Eq. ([IJ in transition metal com- 
plexes has been first discussed in [TT]. Here we reformu- 
late this strategy in terms of projection operators. 

We label Ho the unperturbed Hamiltonian in the ab- 
sence of the time-odd fields B and fi 1 . The well-known 
microscopic operators |12j . linear and bilinear in B and 
[i 1 , describing the interaction between electrons and 
fields are given by [13] (i) the Zeeman term J^ z = — fi e ■ B 
where \i e is the electronic magnetic moment and (ii) the 
four Ramsey terms, describing the interaction between 
the electrons and fi 1 . The Ramsey terms can be further 
partitioned in a contribution linear in fj, 1 , and a contri- 
bution that is bilinear in B and fj, 1 . We collect the terms 
linear in fi 1 in J-^p,^, where the time-odd operator J-J 3 — 
J- 1 p + J~2 fi +-^3 /3 i s t ne sum °f the nuclear-orbit (T[ /3 ), 
Fermi-contact (F^p) and spin-dipolar (F^p) terms [15] . 
The bilinear term is B a T> a p pA, where the time-even op- 
erator V^a describes diamagnetic shielding in closed-shell 
molecules [51 IT5]. 

The thermally averaged field at the position of nucleus 
I is obtained by averaging the total hyperfine field: 
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where p — exp[— /3(Hrj +Jf%s)] is the density operator and 
f3 = (k-QT)^ 1 . Next, we expand p to first order in 
QTJ [T3] and retain those terms of Eq. ^ that are linear 
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in B a . The result is expressed as a sum over the energy 
levels n (with degeneracy w n ) of Hq: 
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where P„ = Y^i=i \ni)(ni\ is the projector on level n, 
Q„ = 1 — P n , and Z = X) n ^ne~' 3e ". This general ex- 
pression for the shielding tensor formally resembles the 
Van Vleck equation for magnetic susceptibility |15) . 

A common situation consists of a system with only 
one thermally occupied degenerate energy level. The 
shielding tensor of such system can be decomposed as 
(er) = o- p +<7 r , where the index / has been dropped to ease 
the notation. The temperature-independent term a 1 is 
universal to degenerate and non-degenerate states alike, 
and is known as the Ramsey shielding in the latter case. 
The effect of degeneracy is to shift cr r by a temperature- 
dependent amount o~ p , the paramagnetic shift, which will 
be the focus of this Letter, and which, from Eq. ([3]), can 
be expressed as (Pq is the projector on the ground state): 
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In Ref. [7J, 

a a/3 was calculated from the parameters of 
an EPR spin Hamiltonian, expressed in the usual way 
as /ieS • g • B + S • A • I. However, this approach is 
only correct if 5 < 1, i.e., for 3- fold and 2- fold electronic 
degeneracies, but it is only an approximation for higher 
spin. [16] Moreover, while the approximation works for 
pure spin states, definition of spin Hamiltonians in the 
strong SOC limit needs additional caution. 

To describe the strong coupling case we have in fact 
to resort to the concept of fictitious spin, and regard 
the degenerate manifold as a spin multiplet 5, so that 
25 + 1 = uj. Once the ui wavefunctions of the ground 
manifold have been optimized via ab initio methods, the 
lu x uj matrix representation X of any operator X in the 
basis of these wavefunctions can be reproduced by an 
effective operator in the spin space [16j . sum over irre- 

(k) 

ducible tensor operators S q of rank k, with spherical 
components q — —k, . . . , k labelled by Latin indices: 
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Here the matrices S^, representation of 5_^ on the fic- 
titious spin basis, form an orthogonal basis for the vector 
space of all complex square matrices of dimension 25-1-1. 

(k) 

The coefficients Xq can thus be determined from the ab 
initio matrices X by orthogonal projection: 
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Our first aim is to reformulate the Curie shielding ten- 
sor Eq. Q in terms of spin Hamiltonian parameters valid 
in the strong SOC limit, as these parameters can be rou- 
tinely measured in EPR, and computed via accurate ab 
initio methods. To this end, we apply Eq. Q to spin- 
decompose the three components of the microscopic Zee- 
man and hyperfine "fields" appearing in Eq. Q, leading 
to three sets of X-numbers for each field. A numerical 
example of how a projection is carried out using Eq. ^ 
and Eq. ^ is reported in the EPAPS [25] . 

Following usual notation we name the resulting gen- 
eralized EPR "tensors" collecting these sets of numbers 
g q k a (Zeeman) and A q k 2 (hyperfine). For k — 1, we obtain 
the usual EPR g-tensor and A-tensor previously consid- 
ered OEJ. For a general degeneracy, and arbitrary SOC 
strength, we obtain: 

P He, a P0 = -MB E(-!) 9 ^ 42 ( 7 ) 
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with k < 2S and odd [TT], and the ^-factor of the 
nucleus in question. 

Substituting these expressions in Eq. Q yields 
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This expression generalizes all formulas previously pro- 
posed in the literature to the strong spin-orbit coupling 
limit, and to arbitrary degeneracy, thus representing one 
of the main results of this Letter. Interestingly, only 
products of same-rank tensors enter Eq. Q. 

Next, we look for an expression for cr p that is more 
general than Eq. |[9j), accounting for a time-even pertur- 
bation Hi that weakly splits the 25+1 degeneracy, which 
represents a common situation (see Figure 1). In cubic 
symmetry Hi may describe the vibronic coupling of a 
Tg electronic state to a Jahn- Teller active mode. For a 
pure spin degeneracy (5 > 1) [7j, Hi leads to the EPR 
zero-field splitting (ZFS) Hamiltonian. 

The same arguments [16] leading to Eq. |7]) and Eq. ^ 
can be used to project Hi on the ground-manifold, 
irrespective of its microscopic origin, and map it 
into a spin Hamiltonian. This leads to PoHi-Po = 
J2k q{~ lyDfF 1 S_q, where k is an even integer [UJ (k < 



(k) 

25), and the generalized "ZFS tensors" D q ' are eval- 
uated via Eq. For instance, the projector Pq on 
an ab initio r§ ground manifold is Pq = J2q\l)(l\^ 
where q — K, X, jj,, v [18]. The rotational properties 
of a spin basis quantized along a C4 axis (Z) readily 
lead to the mapping [T5]: \k) -> |3/2), |A) -> |l/2), 
— > | — 1/2), \v) — > I —3/2) , and, for an axial distortion, 
to the spin Hamiltonian PqHiPq = D{S\ — 5/4), where 
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FIG. 1. Scheme of an axially split Tg ground state and as- 
sociated projectors (see text). Left: Fg ground state (GS) of 
Ho, with ab initio projector Po mapped to a spin (S = 3/2) 
projector Pg = 3/2- Right: splitting induced by Hi, described 
in GS by the spin Hamiltonian PoHiPo = DS%, whose eigen- 
functions define the split projectors 0±3/2 = 3/2) (3/2 + | — 
3/2){-3/2|,e±i /2 = 1 1/2) (1/2| + | -1/2) (-1/2| . 



D = (k|Hi|k) = (v|Hi|i/) = — (A|Hi|A) = -(m|Hi|/x). 
The paramagnetic shift for a split manifold can now be 
obtained only using the spin basis, via three simple steps. 

First, we diagonalize the u> x uj spin Hamiltonian 
PoHiPo, and find the states \ip\ a ) an d split energies e\, 
where A = 1, . . . , u/ (u/ < 25 + 1) labels the new low- 
lying degenerate manifolds, while a spans a single man- 
ifold. Next, the eigenstates \ij}\ a ) are used to build new 
spin projectors a = J2 a lV-'Aa)(V'Aa|, each associated to 
a degenerate manifold A. For the axially split r§ case, 
as illustrated in Figure 1, DS\ is diagonal, T$ splits into 
2 Kramers doublets (KD), and the new spin projectors 
can be written as (D < 0) 6±3/2 for the ground KD, 
and 0±i/2 for the excited KD (see caption to Figure 
1). Finally, the third step to obtain the paramagnetic 
shift for the split case is to apply Eq. but instead 
of microscopic Hamiltonians and projectors, we now use 
the generalized Zeeman and hyperfine spin Hamiltonians 
{PoHe,aPo an d PsJ-pPo), and, crucially, the spin projec- 
tors defined on the split manifold (0a)- 

Following this strategy the paramagnetic shift of a split 
manifold can be solely expressed as function of g( k > and 
tensors, the rank k referring to the (25+l)-manifold: 
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where the partition function is defined over the split man- 
ifold as Zo = 2^a w Ae~^ £A , and the non-trivial tempera- 
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Note that Eq. ( 10 ) is different from the ZFS-expression 
proposed in Ref7[7 , in that (i) it depends on all possi- 
ble products of tensors of all ranks allowed by the size of 



the degenerate manifold (ii) the "Curie term", Eq. (Ill, 
cannot be written in terms of averaged products of spin 



operators as (Sq k ^£& (iii) Eq. (10) contains a new "or- 



V 



bital" term, Eq. (12), giving rise to a non-Curie temper- 



ature dependence. The new term is essential to achieve 
the correct high-temperature limit, i.e. Eq. Eq. (10), 



generalization of Eq. ^ to weak-split degeneracies, is the 
other main result of this Letter. 

Applications. We illustrate the theory by evaluating 
the NMR paramagnetic shift of strong spin-orbit coupled 
dopant lanthanide nuclei in fluorite crystals, generally de- 
noted as CaF 2 :Ln. Among the Ln-ions occupying a cubic 
crystal site Pr 2+ , Nd 3+ , Sm 3+ , and Dy 3+ are known to 
have a 4- fold degenerate ground state [HO [25] , with sym- 
metry T% in the cubic double group O* . Our theory is 
essential to correctly describe NMR for these nuclei. 

Associating the r 8 states to a spin 5 = 3/2 as seen 
previously, from Eq. Q and Eq. ([8| one obtains: 
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are the components of that transform as a vector 
under octahedral symmetry |20j . Symmetry reduces the 
number of non-zero tensor components to just four, de- 
noted by g, g', A, and A'. Next, by means of Eq. the 
paramagnetic shielding tensor (diagonal and isotropic in 
this symmetry) can be expressed as 
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To calculate <r p via Eq. (16), we need the relevant 
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rank-1 (g and A) and rank-3 (g' and A') tensor com- 
ponents. We choose the multiconfigurational wave func- 
tion method CASSCF, followed by a non-perturbative 
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FIG. 2. (color online) Plot of a p T vs. T for the Nd impurity 
nucleus in CaF2:Nd crystal computed using the values in 
Table [I] The solid (blue) line corresponds to a p T for an exact 
F$ ground state. The dotted (red) and dashed (green) curves 
correspond to a p T and er^T, respectively, for an axially split 

(D = — 5cm -1 ) Fg state, calculated using Eq. ( 17 1. 



calculation of the SOC as implemented in the code MOL- 
CAS |21j . The Ln- impurity and its first coordination 
sphere (LnFg - ) are treated as an ab initio fragment em- 
bedded in the CaF2 crystal potential, which is described 
by a combination of ab initio model potentials (AIMP) 
[22) and point charges in a finite cube centered on the 
Ln-ion (see EPAPS [23] for further details). 

The ab initio matrix elements of the magnetic moment 
in the T% ground state were used to determine g |24j and 
g' . Unfortunately, the calculation of the hyperfine con- 
stants is currently out of reach of the CASSCF method 
for most systems. Hence we chose to evaluate A and A' 
by assuming that A/g = A' / g' = Aj/gj, where gj and 
Aj are the free-ion Lande factor and hyperfine constant. 
This relation is exact in the limit of negligible J-mixing 
|25j . arguably a good approximation for Ln-complexes. 
Estimating Aj/gj from tabulated experiments [231 121)] 
allows us to obtain A and A' from the computed values 
of g and g' . All results are in Table [I] The ab initio g 
values for Pr 2+ and Dy 3+ compare reasonably well with 
available EPR data: g cxp = -0.881 and g' exp = -0.173 
(Pr 2 +) [28, «? exp = 4.4 and <4 p = 0.19 (Dy 3 +) JTS]. Note 
that contributions arising from higher order EPR tensors 
(g' and A') in CaF 2 :Nd 3+ (CaF 2 :Pr 2 +) account for 54% 
(46%) of the shielding tensor [27] . 

Finally, if the cubic environment of the impurity un- 
dergoes Jahn- Teller axial distortion along z, described by 
the ZFS fictitious spin Hamiltonian D(S^ — 5/4), appli- 
cation of Eq. (10 1, leads to the following components of 



TABLE I. Ab initio (CASSCF) EPR parameters and NMR 
paramagnetic shift in the Fs ground state of CaF2:Ln n+ , com- 
pared with available experimental data. The hyperfine con- 
stants (10 -2 cm -1 ) are divided by the nuclear (?n factor to 
make them isotope-independent. Values of Ta p (K) are ob- 
tained from Eq. (161. In parenthesis we report the ab initio 



shift that would be obtained if g' and A 1 were neglected. 





g 


g 


A/gN 


A'/gs 


Ta p 


Tcr p (exp) 


p r 2+ 


-0.799 


-0.245 


-2.17 


-0.666 


-106 (-57) 


-94^] 


Nd 3+ 


-0.802 


-0.293 


-2.66 


-0.971 


-155 (-71) 




Sm 3 + 


-0.233 


-0.0632 


-2.77 


-0.754 


-36 (-21) 




Dy 3 + 


4.257 


0.198 


6.05 


0.281 


-867 (-850) 


-90fQ 



Calculated from EPR data reported in Ref. 19 26 



the anisotropic Ln-shielding tensor in the split-Fg 



state: 
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Note that in Eq. (17) products of tensors of different 
ranks appear, and the last term of a\ is the new "orbital" 
contribution Eq. (12 1. Figure 2 shows the non-trivial 
temperature dependence of Cy and for the weakly 
split case, compared to the simple Curie term for a Tg 
ground state, in CaF 2 :Nd 3+ . This non-Curie behavior 
provides an experimental strategy to probe small inter- 
actions between nanomagnet and environment, by moni- 
toring the paramagnetic shift as function of temperature. 

In conclusion, we presented a general theory of NMR 
paramagnetic shifts, valid for arbitrarily degenerate elec- 
tronic states. The shielding is a function of general- 
ized EPR Zeeman and hyperfine tensors of all symmetry- 
allowed ranks. The significance of the approach has been 
demonstrated by high-level ab initio calculations. We be- 
lieve the theory will be valuable to interpret pNMR data 
for lanthanide-based magnetic materials, and for the de- 
velopment of ab initio methods for pNMR. 
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